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H-I Abstract 
OJQ, 

We give a perturbative analysis for an infinitesimally thin cylindrical shell com- 
' posed of counter rotating collisionless particles, originally devised by Apostolatos and 

\ Thorne. They found a static solution of the shell and concluded by C-energy argu- 

^ . ment that it is stable. Recently, the present authors and Ida reanalyzed this system 

by evaluating the C-energy on the future null infinity and found that the system has 
an instability, though it was not shown how the system is unstable. In this paper, it 
is shown in the framework of the linear perturbation theory that, if the constituent 
particles move slowly, the static shell is unstable in the sense that the perturbation of 
its circumferential radius oscillates with exponentially growing amplitude, whereas if 
the speed of the constituent particle exceeds a critical value, the shell just expands or 
contracts exponentially with time. 
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1. Introduction 



Nowadays, it may be a matter of time before one makes direct observations of gravi- 
tational waves, which will be a great success in experimental physics. Thus, it becomes 
more important to theoretically understand well about spacetime dynamics accompanied 
with emission of the gravitational radiation. However, it is not so easy task, since it is, in 
general, difficult to solve the Einstein equations because of their nonlinearity. One efficient 
way is to impose a symmetry on spacetime. The simplest one is the spherical symmetry, but 
there is no freedom of gravitational radiation in spherically symmetric system, and hence 
this assumption is not adequate to study the effects of gravitational radiation. As another 
simple symmetry, cylindrical symmetry has been sometimes considered. The cylindrically 
symmetric spacetime has a degree of freedom of gravitational radiation known as Einstein- 
Rosen gravitational waves] The cylindrical system has been studied in connection with 
the gravitational waves by some authorsP'®'^ 

Apostolatos and Thorne studied an infinitesimally thin cylindrical shell composed of 
counter rotating collisionless particles with vanishing total angular momentumP Hereafter, 
we refer to this as the Apostolatos- Thorne(AT)-shell. The AT-shell admits a static config- 
uration, and Apostolatos and Thorne investigated its stability by the argument based on 
the C-energy which is a quasi-local energy per unit length in the translationally invariant 
direction.^ They concluded that the static configuration of the AT-shell is stable, and any 
dynamical configuration will finally settle down into the static state by releasing the energy 
through the gravitational emission. By contrast, the present authors and IdsP have shown 
that some class of momentarily static and radiation free initial configurations does not settle 
down into static equilibrium configurations, or otherwise infinite amount of the C-energy 
through the gravitational radiation is released to the future null infinity. It implies the ex- 
istence of an instability in this system. It should be noted that the C-energy argument by 
Apostolatos and Thorne is restricted in a bounded domain, whereas the argument in Ref. [8] 
is based on the conservation of the C-energy in the infinite domain. 

Recently, Gleiser and Ramirez® have studied some aspect of the relativistic dynamics 
of the AT-shell. They have analytically solved the equation of motion for the shell and the 
Einstein equations by imposing a condition that the interior of the AT-shell is always flat. 
They have also investigated a perturbation around the static configuration and found that 
there are oscillating solutions for any chosen period. It should be noted that they analyzed 
the cases which permit the C-energy input to the shell from the past null infinity and is 
different from the situation considered in Ref. 

In this paper, we reanalyze the AT-shell model in the framework of the linear perturbation 
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theory. We are interested in stability of the AT-shell in the case that there is no C-energy 
input to the shell from infinity, and hence we assume that there is no incoming gravitational 
radiation from past null infinity. 

Hamity, Cecere and BarraccP^ discussed the stability of the static state of the shell, and 
concluded that there are both stable and unstable static configurations for the AT-shell. We 
will mention their result in connection with our present study in the last section. 

This paper is organized as follows. In section El we describe the geometry of the cylindrical 
shell and construct the AT-shell model by using the coordinate system comoving to the AT- 
shell. In section [31 we describe the static solution for the shell. In section HI we study the 
stability of the static AT shell by a perturbative analysis. Finally, a brief summary of the 
main results is given in section [5j 

We adopt the geometrized units c = G = 1. 

§2. Cylindrical shell composed of counter-rotating particles 

In this section, we give a description of the AT-shell model which is, as mentioned, an 
infinitesimally thin cylindrical shell composed of counter rotating collisionless particlesP 
For our present purpose, we adopt the coordinate system comoving to the shell, in which the 
spatial coordinates of the shell are kept constant. This coordinate system is different from 
the original one given by Apostolatos and ThorneP 



t 




Fig. 1 . The world volume U of the AT-shell divides the spacetime into two regions M± . 

The AT-shell divides the spacetime into two regions: one is the inside of the AT-shell 
M~ and the other is the outside of it M + (see Fig. [T]). By denoting the world volume of the 
AT-shell by E, the whole spacetime is given by M = M~ U E U M + . Both regions M ± have 
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the whole cylinder symmetry, or in other words, the line element is written in the form, 

ds 2 = e 2 ^\-dt 2 + dr 2 ) + f3 2 e- 2 ^d<t> 2 + e^dz 2 , (2-1) 



where 7, ip and (3 are functions that depend only on t and r. We assume that both regions 
M ± are vacuum. The domains of coordinates are — 00 < t < 00, — oo<2;<oo,0<r<oo, 
and < < 2ir. The axis of rotational symmetry is located at r = 0. The coordinate bases 
d/d(f) and d/dz are the rotational and translational Killing vectors, respectively. 

The Einstein equations in the vacuum region lead to the equation for /3, ip and 7 as 



P"-P = 0, (2-2) 

W-O + /^-/3ty' = 0, ( 2 ' 3 ) 

(f3 2 - (3' 2 )i + (3(3' \^' 2 + - 2P$W>' + iP'P" ~ PP 1 ) = (24) 

((3' 2 - /3 2 ) 7 + (3^ 12 + ^ 2 ) - 2^'#' + ((3(3" - f3'/3') = (2-5) 

7 - 7" + - ^ /2 ) = 0, (2-6) 



where A and A' mean derivative of A with respect to t and r, respectively. In the original 
description of the AT shell, (3 is set to be r, i.e., (3 = r. By contrast, we do not set it here 
so as to introduce the coordinate system comoving to the shell later. 

The stress-energy tensor of the AT-shell is infinite since the finite energy is confined 
within the infinitesimally thin region. The Ricci tensor diverges on the AT-shell through 
the Einstein equations, and hence the AT shell is the so-called s.p. curvature singularity.^ 
However, since the co-dimension of a shell is one, this singularity is so weak that the metric 
is defined on the AT-shell. Hence, even if the spacetime is singular on the AT-shell, it 
can consistently be treated by the Israel's metric junction met ho d.^'^'^ The junction 
conditions between M ~ and M + at E require the continuity of the metric and specify the 
discontinuity of the extrinsic curvature K^p compatible with the stress-energy tensor of the 
AT-shell. 

We require that the Killing vectors d/d<p and d/dz are continuous at E, and this re- 
quirement leads to the continuities of the metric functions (3 and ip. We do not require the 
continuity of 7, and this means that the coordinate function t and accordingly, the coordinate 
basis, d/dt, may not be continuous at the E. Although we can require that the coordinate 
function t and accordingly, the coordinate basis, d/dt is also continuous at E, we will not 
do so in this paper. This means that 7 may not be continuous at the AT-shell. The reason 
of this choice is that the solutions for the metric functions do not take simple forms in the 
continuous time coordinate, whereas the discontinuous one makes those solutions simple. 
The time coordinate for the interior region is denoted by t_, whereas that for the exterior 
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region is denoted by t+. The radial coordinate at the AT-shell is denoted by R(t), where r 
is the proper time of an observer at rest on the shell. 

We introduce the proper reference frame of an observer riding on the AT-shell as follows: 

:our velocity of the shell, (2-7) 

outward unit vector normal to the shell, (2-8) 

(2-9) 

(2-10) 







e (T) = 


d -x 9 +v d 

dr dt± dr 






e ( „) = 


d -v 9 hi 9 

dn dt± dr 




1 d 


e «>) = 


fie-* dtf 




1 d 


e W = 





where 



dR . dt- 



V:=^P, and X ± := ^ = vV^-VO + V\ (2-11) 
dr dr 

The subscripts or superscripts + and — are used to denote quantities evaluated on the outer 
and inner faces of the shell, respectively. 

Each constituent particle of the AT-shell has an identical rest mass. Half of the particles 
orbit around the symmetry axis in a right-handed direction with angular momentum per 
unit mass a, and the other half orbit in the opposite, left-handed direction with angular 
momentum per unit rest mass —a. Therefore, the net angular momentum of the AT-shell is 
zero. The absolute value of the specific linear momentum u of a constituent particle, which 
is defined by 

is an important quantity to describe the AT-shell. The other quantity characterizing the AT- 
shell is the shell's rest mass per unit translational Killing length, which is denoted by A. The 
rest mass of each constituent particle is conserved, so is A. Then, the surface stress-energy 
tensor of the AT-shell is given by 

S = cre (T ) <g> e (r ) + Te w <g> e w , (2-13) 

where we have defined the energy per unit area and the surface stress of the AT-shell as 

a:= XVTTV \ and T := ^ , (2-14) 

respectively. 
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Israel has shown that the Einstein equations on the AT-shell reduce to 



where K aj3 is the extrinsic curvature of the shell's outer and inner face, and h a p is the induced 
metric on the E. 

Now, we impose the comoving condition for the AT-shell, i.e., V = 0. The comoving 
condition fixes the metric function /3, but it does not necessarily imply (3 = r. In this coor- 
dinate system, the coordinate radius R at the AT-shell does not characterize the evolution 
of the shell, but its circumferential radius 71 defined by 

U := Pe-%= R , (246) 

does instead of R. 

In the comoving coordinate, the non- vanishing components of the extrinsic curvature are 
given by 

%) = ^n-7,n> (2-17) 

K {m) = (ln/3), n -^ n , (2-18) 

#(,)(*) = tn- (2-19) 

Therefore, the Israel's junction conditions give the following equations: 

[.-VV-,f^^, (2.20, 



e" V ( | - if/ 



2A(1 + 2u 2 ) 



(2-21) 



[e- ? eV]S- / ^_ , (2-22) 



/3vT 



where 

[A]+ = A + -A". (2-23) 
Using the continuity of ip and /3 at Z", we can rewrite the above equations as 

^V]! = -^=e^, (2-24) 



PV1 



u 



[e-V] + = ^^e^, (2-25) 



[ e -f/3'] _ = -4Av / IT^e- ,/ '. (2-26) 
These are the junction conditions under the comoving condition. 
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§3. static solution 



As was found by Apostolatos and Thorne,® there is a static solution for the AT-shell 
which is regular except at the shell itself. In the interior region, the metric functions of the 
static solution are 

P~ = r, V~ = V>., 7" = 0, (3-1) 
where ip s is a constant. In the exterior region, the metric functions are 

/3 + = r, ifj+ = ijj s -2y\n^ 7 + = 2 ln(l + 2y) + Ay 2 In ^, (3-2) 
where we have defined a parameter 

y ■= ( 3 - 3 ) 



The junction conditions require the relation 



Ae"* = (3-4) 
(l + 2y) 2 - [6V 

Since we may set xjj s = by rescaling the coordinates t, r and z, the static solutions form a 
one parameter family. If we fix y, then the solution is uniquely determined, and hence y is 
a good parameter to characterize the static solutions. 

§4. Perturbative analysis around the static solution 

In this section, we give a perturbative analysis of the static AT-shell. If there is a smooth 
one-parameter family of dynamical exact solutions for the AT-shell, which includes a static 
solution as its member, we can construct dynamical solutions from the static solution by the 
perturbative method.^ Here, we assume that such a one-parameter family exists. Denoting 
the parameter by e, we write the metric functions in the form 

/3(t,r;e) = ^ n /3 n (t,r), (4-1) 

n=0 

^(t,r;e) = ^5> n (t,r), (4-2) 

n=0 

7 (t,r;£) = 5>" Tn (t,r). (4-3) 



n=0 



The static solution is obtained by putting e = 0, or in other words, (3 , ip Q and 7 are 
functions given by ( 13 - 1 p and f)3-2p for the interior and exterior regions, respectively. We can 
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derive the Einstein equations of each order with respect to e, and solve them by successive 
approximation up to the order we need. In this paper, we solve the equations up to 0(e), 



i.e., the linearized Einstein equations in the static background. 
The Einstein equations of 0(e) are 

Pi ~ h = 0, (4-4) 

M'i + - Mi + M% + PW = o, (4-5) 

Ti " 7i + 2^>i = 0, (4-6) 

p' Q ft - 2(3'm + m' 2 p[ + pi& x + p , Q ^p 1 - /# 7 ; + 2^^^ = 0, (4-7) 

P'oPi - Mo Pi - P'qPi - (3' 2 ii + 2PoPMi = 0. (4-8) 



Eqs. ()4-4|) - (j4-6p are the evolution equations for Pi, ipi and 71, whereas Eqs. (14- 7p and (14- 8p 
are the constraint equations which include no second order time derivative. 

In order to study the stability of the static AT-shell, we assume that there is no energy 
input into the AT-shell. Therefore we impose the outgoing wave boundary condition in 
which there is no ingoing gravitational wave from past null infinity. In the interior region, 
we impose the regularity condition at the symmetric axis. Firstly, we solve these equations 
in the interior and exterior regions by imposing the regularity condition at the symmetry 
axis and the outgoing wave boundary condition, and then impose the junction conditions to 
glue these solutions at E. 

4.1. Solutions for the interior region 

By using the expression of the static solutions in the interior region (13-1 1) , we rewrite the 
linearized Einstein equations for the interior region as 



pf - pi = 


(4-9) 


i^i" + ~i>i - $T = °> 

r 


(4-10) 


li" ~ li = 0, 


(4-11) 


a-" -' 
Pi =7i , 


(4-12) 


k' = ii- 


(4-13) 



The last two equations can be integrated to give the relation 

P~' + C- = ^, (4-14) 

where C~ is an integration constant. Therefore, 7^ is easily obtained once is known. 
The others are three homogeneous equations. Since the background solution is static, we 
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can assume that all perturbations have a time dependence e~ luJ ~ l ~ and solve for U- in order 
to know whether complex or pure imaginary frequencies exist. The solutions for the three 
homogeneous equations are given by 



-iuj—t- 



A? 



tu_)e iuJ ~ r + B p ~ (w_)e 
(i) 



A^ {uJ)H^'{uj-r 



= Re 
ifti = Re 
lx = Re 

where A p ~ , B p ~ , A^' and B p ~ are functions of 

(2) 

Hq are the Hankel functions of the first and second kind of order 0, respectively. 

Now, we impose the regularity condition at the symmetric axis, r = 0, in accordance with 
HaywardP^ 1 The regularity condition can be expressed in a geometrically invariant manner 
by using the norm of the Killing vectors, 



(4-15) 
(4-16) 
(4-17) 

'Re" denotes the real part, and Hq\ 



p:=\d4,\=Pe-*, £:=\d z 



(4-18) 



At the symmetric axis, p should vanish. The axis is said to be regular if £ is finite and 



Vp ■ Vp = 1 + 0(p 2 ) 
Vp-W = 0(p), 
W- W = £>(1), 



(4-19) 
(4-20) 
(4-21) 



where V is the covariant derivative of the spacetime metric. These regularity conditions are 
satisfied if and only if 



A 13 +B 13 = 0, A^ - = 0, C 
Therefore, we obtain the regular interior solutions as 



0. 



ft- 


= Re 






= Re 


2e~ iuJ 


7i" 


= Re 





t -A?-{u:J)( 

-*-^"(w_)Jo(a;_r) 
A p ~ (e iuJ ~ r + 



) 



(4-22) 



(4-23) 
(4-24) 
(4-25) 



where Jo is the Bessel function of order 0. If uj- is real and negative, we have to replace 
Jo(u^r) by Jo(|o;_|r). However, since we are interested in complex or pure imaginary solu- 
tions for oj_, we need not treat such a case. 
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4.2. Solutions for the exterior region 

Using the fact that ip'^ = — 7q, the constraint equations (14 -7p and (14- 8[) become 

Pt" - {it dt)' + 2/3 + ^o + >i + ' - /?o + V = 0, 
Pt' ~ it'Pt + - Pi' it = o. 



(4-26) 
(4-27) 



By noting that /3qiPq is constant, we can integrate the above two equations and obtain 



it = -^7 ( Pi - it'Pt + W^t'rt ) + 



(4-28) 



where C + is an integration constant. By using the background solutions (I3-2|) . the evolution 
equations for 0i, ipi and 7^ become 



at" - ft = 0, 

jf - 7+ + 2^+Vf = 0. 



(4-29) 
(4-30) 
(4-31) 



Though the equation for ipf is an inhomogeneous equation, one can easily see that tpi fit is 
a particular solution to the equation^]. By the same reason as the interior solutions, we can 
assume the time dependence e~ lu)+t+ for the exterior solutions. The outgoing wave solutions 
for the above equations are written as 



^i 


= Re j 




4+ 




2yPt 




J 

r 


it 


= Pi' ■+ 





iu + (t + -r) 



(4-32) 
(4-33) 
(4-34) 



where we have introduced the outgoing homogeneous solution for Eq. (14-301) as 

g} ■= ^ 



(4-35) 



If uj + takes a negative real value, we have to replace H^iou+r) by H^\\u+\r) so as to keep 
the outgoing wave condition. However, since we are interested in complex or pure imaginary 
solutions for oj + , we need not treat such a case. 



Note that ip^ = —2y/r from Eq. 



10 



4.3. Solving the junction conditions at E 

Now, we glue the interior solution to the exterior one at E. For that purpose, we have 
to relate the time coordinates t_ and t+ at E. Firstly, we relate t- and t + to r. From 
Eq. ( 12- lip and the comoving condition V = 0, we have 



~dr 



-7 +i ) 



+ ee^(^-^) + 0(e 2 ). 



(4-36) 



The second term in the right hand side in the above equation contributes in all relevant 
equations only of 0(e 2 ), since the background solution does not depend on Therefore, 
we may set, in the linearized equations, t_ = e^ s r. Similarly, for the exterior time coordinate, 
we may set t + — (1 + 2y)~ 2 e^' s r. From these relations, we obtain 



*+ = (! + 2?/)~ 2 t_ on E. 



(4-37) 



The continuity of /3 and ip at E gives 



01 \r=R 
i>T\r=R 



ft | r =R, 
1pf\r=R- 



(4-38) 
(4-39) 



The above conditions imply that the time dependences of solutions for ft and ip\ in M± 
are the same at E, and hence ou + t + = should be satisfied. Then from Eq. (j4-37j) . 

we have w + = (1 + 2y) 2 u_. Hereafter, for notational simplicity, we abbreviate w_ as oj, or 
equivalently, 



0J- 



u+ = (l + 2y) 2 u. 



(4-40) 
(4-41) 



The junction conditions of 0(e) are obtained by directly expanding Eqs. (I2-24|) . (I2 25[) 
and (12-261) with respect to e as 



[e- 70 (7;-7i7^)]! 



■2y 



1 ft~ , 1 + 22/ 



R[l + 2yf \l + y R l + y 
Ay 2 



1 3 + 22/ft" 



i 



4y(l + y) 
(l + 2j/)2 Vl + y 



1 + 2/ i? 

_y_P\_ 

l + y R 



r=R 



r=R 



r=R 



(4-42) 

(4-43) 
(4-44) 



Since ft and are continuous at r = i?, one may replace ft and in the right hand 
side of the above equations with ft 1 " and t/^, respectively. By substituting Eq. (14- 14[) with 
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C~ = and Eq. fl4-28j) into Eq. fl4-44j> . and by using the static solutions (13d]) and (J3Z2J), we 
obtain 

C + = 0. (445) 
Then, four equations (!4-38j) . (!4-39p . ()4-42p and ( I4-43P determine four amplitudes and 

The continuity conditions (I4-38P and (I4-39P give 

- AP + A 13 ' 

(l + 2y) 2 e tU}R —— -2i——sinuR=0, (4-46) 
R R 

- A^ 

(1 + 2y) 2 A^ + H^\ujR) - 2y(l + 2y) 2 e iuR — - 2A^ J Q (uR) = 0, (4-47) 

R 

where, for notational simplicity, we have introduced a new quantity R — (1 + 2y) 2 R. The 
junction conditions ( l4-42p and (I4-43P lead 

-uR(l + 2y) A A" p+ H^ ) {LoR) - 8y 2 (l + 2y) 2 A^ + H^\coR) + 2(1 + 2y) 2 uRA^ J^loR) 
+2y(l + 4y 2 )(l + 2yfe^ R ^- - sin coR - ^ + ^ ^ J M) = 0, 

(4-48) 



and 



- 4' 3+ - 
-w 2 i2 2 (l + 2?/)V^— - - 4?/ 2 (l + Ay 2 )(l + 2y) 2 e luR — + 4y(l + 2y) A uRA^ ' h[ 1] } (uR) 

R R 

+lQy 3 (l + 2y) 2 A^ + H^\uR) + 2i(l + 2y) Vtf 2 ^- sinwi? - J^—A^J^R) 

R (1 + 2/ J 

8n/ 2 (3 + 2y) A/ 3 " „ ,„ „ . 

+ ( H») T mil ° c ' (449) 

respectively. These equations can be put into the following matrix form: 

4 

M a Aj = o, (4- 5 °) 

where z = 1,2,3,4 and we have defined 

A X :=A^\ A 2 :=—, A 3 := A^~ , A 4 := — — . (4-51) 
The components of the matrix My are given by 

M n = M13 = M 24 = 0, (4-52) 
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Mi 2 = (1 + 2y) 2 e iujR , (4-53) 

M 14 = -2isinuR, (4-54) 

M 21 = (l + 2y) 2 if ( S 1) (a; J R), (4-55) 

M 22 = -2y(l + 2y) 2 e™ R , (4-56) 

M 23 = -2J {uR), (4-57) 

M 31 = + 2y) 4 H[ 1 \ujR) - 8y 2 (l + 2y) 2 H^\ujR), (4-58) 

M 32 = 23/(1 + 4y 2 )(l + 2yfe wk , (4-59) 

M 33 = 2(1 + 2y) 2 u; J RJ 1 M?) - ^±M j ( u R) t (4-60) 

(1 + 2/) 

M 34 = - %V , sinwi?, (4-61) 
(1 + 2/) 

M 41 = 4y(l + 2yfujRH i i\ujR) + 16y 3 (l + 2y) 2 H { ^\uR), (4-62) 

M 42 = -w 2 i? 2 (l + 2y) 6 e^ - 4y 2 (l + 4y 2 )(l + 2y) 2 e i " R , (4-63) 

M 43 = -J*f—J ( U R), (4-64) 
(1 + 2/) 

82 U 2 ('3 + 2v) 

M 44 = 2i(l + 2y) 2 u 2 R 2 sinuR + ^ + ^ sinwA (4-65) 

In order that there are non-trivial solutions for Eq. ( l4-5Up . the following equation should be 
satisfied; 

det (Mij) = 0. (4-66) 

For fixed y and R, this is an equation for frequency u. Here, instead of f)4-66p . we consider 
a real equation 

| det (M i:j ) | 2 = 0. (4-67) 

Using a dimensionless quantity z = uR, Eq. (14-671) is rewritten in the form 

256y 2 (l + 2y) 12 



where F(z) is defined by 



-F(z) sin 2 z = 0, (4-68) 



F(z) 



2y{\ + y) ((1 + Qy + Ay 2 )z 2 - 2y) J (z) 
+z (y(2 + y(l + 2y) 2 ) - (1 + y) 2 (l + 2y) 2 z 2 ) J^z)] H® ((1 + 2y) 2 z) 
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+ (1 + 2y) 2 z[{z 2 + y(2 + y)(z 2 - l))J Q (z) + 2y(l + y)zJ l (z)jH{ 1) ((1 + 2y) 2 z) 

2y{\ + y) ((1 + Qy + 4y 2 )z 2 - 2j/) J (z) 

+z (y(2 + y(l + 2y) 2 ) - (1 + y) 2 (l + 2y) V) J^)) tf< 2) ((1 + 2y) 2 z) 

+ (1 + 2y) 2 z(V + y(2 + y)(z 2 - 1) J (z) + 2j/(l + y)^^)))^ ((1 + 2y) 2 z) 

(4-69) 

It is easily seen that, for integer n, there are infinitely large numbers of real solutions 
z = nil, or equivalently, 

- = (4-70) 

These are radial oscillation modes of the shell around the static configuration. The existence 
of these oscillation modes is compatible with the C-energy argument given by Apostolatos 
and ThorneP The gravitational emission may cause the damping of the oscillation, but this 
effect appears in 0(e 2 ). 

4.4. Unstable modes 

In this subsection, we look for zero points of the function F(z) in complex z plane. We 
solve the equation 

F(z) = 0, (4-71) 

by using "FindRoot" program in Mathematica (version 7.0). As mentioned, the parameter y 
characterizes the static solution, and hence we search for the solutions of u as a function of 

y- 

Since we invoked the numerical method to find the solutions, we have investigated a 
limited number of points in the domain < y < oo, not all y. But, as far as we have 
investigated, there exists at least one unstable mode for each y. Therefore, it is reasonable 
to conclude that there are unstable modes for all y > 0. The angular frequency u of the 
unstable solution is written in the form 

u = u R (y) + iu l (y), (4-72) 

where wr is a real function of y, and u\ is also a real function of y but is positive. We also 
found that if oj = w R + ioj\ is a solution, then u = % — iojj is also a solution. 
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We have found sixth classes of the unstable solutions. We depict and u\ of the first 
class as functions of y in Fig. [2J whereas those of the second, third, fourth, fifth and sixth 
classes are dipicted in Figs. [3HZI respectively. 




become pure imaginary. 

In the case of the first class, the real part ojr increases from zero for < y < 0.6, decreases 
for 0.6 < y < 4.61 and then becomes very tiny value e for y > 4.61. The maximum value 
of o/r is about 0.579/? -1 at y = 0.6. The tiny value e depends on initial values in numerical 
investigation, and therefore e will be a numerical error. Hence, it is reasonable to conclude 
that the numerical solutions with the form of u = e ± iui correspond to pure imaginary 
solutions. The imaginary part ui of the first class increases from zero for < y < 4.61 and 
bifurcates at y ~ 4.61: one sequence monotonically increases and the other monotonically 
decreases for 4. 61 <y. At y ~ 4.61, the derivatives of both cur and oj\ with respect to y are 
discontinuous. 

In the case of the second class, the real part increases for y < 1.4 and becomes very 
tiny value for y > 1.4. By the same reason as in the first class, it is reasonable to conclude 
that these solutions are pure imaginary. The imaginary part ui\ decreases for y < 0.8 and 
increases from 0.89-R" 1 for 0.8 < y, approaching R~ l as y -> oo. The imaginary part uj\ of 
the second class is larger than that of the first class. 

In the case of the third class, the real part wr, monotonically increases and approaches 
as y increases. The imaginary part u\ decreases for y < 0.4 and increases from 0.169-R -1 to 
0.236i?~ 1 for 0.4 <y< 0.7. For 0.7 < y, uj\ decreases and approaches as y increases. 
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Fig. 3. The solutions of the second class are depicted. For y > 1.4, the complex modes in this 
class become pure imaginary. 




Fig. 4. The solutions of the third class are depicted. 
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Fig. 5. The solutions of the fourth class are depicted. 

In the case of fourth class, the real part monotonically increases. The imaginary part 
uj\ monotonically decreases and approaches 0. 

The fifth class has monotonically increasing cur for y < 4.61. At y ~ 4.62, cur becomes 
very tiny value, and for y > 4.62 the solution becomes pure imaginary. The imaginary 
part u\ monotonically decreases to 0.075.R -1 for y < 0.8 and monotonically increases for 
0-8 % y 4.61. At y = 4.62, the solution coincides with that of the first class. Therefore, 
this class seems to connect to the pure imaginary modes in the first class for y > 4.62. 

In the case of sixth class, the real part cjr monotonically increases and approaches 0. 
The imaginary part u\ monotonically decreases to 0.016.R -1 for y < 0.44 and increases to 
0.24i? _1 for 0.44 < y < 0.64. For 0.64 < y, uji decreases and approaches 0. 

In order to understand what the existence of the unstable modes implies, we investigate 
the circumferential radius 1Z of the shell. The perturbation of 0(e) for the circumference 
radius of the AT shell is 



\r=R 



Re 



-iujt— 



(e - ^ +y)e iuk A^(uj) 



RH^iu 



(4-73) 



17 



cdR 




i i i i i i i i i i 

2 4 

y 

. 6. The solutions of the fifth class are depicted. At y ~ 4.62, the solution in this class becomes 
pure imaginary and this class seems to connect to the pure imaginary modes in the first class. 
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Fig. 7. The solutions of the sixth class are depicted. 
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From the above equation, we can see that TZi is written symbolically in the form 

Tlx = Re [e~ lult - (d + iC 2 )] = Jcf + C\ e mt ~ cos(u; if t„ - 0), (4-74) 

where C\ and C 2 are real numbers, and we have introduced a constant <fi := tan -1 C 2 jC\. 
The above equation implies that, if cjr does not vanish, the circumference radius oscillates 
and its amplitude grows exponentially with time. If cjr, vanishes, there is a mode of which 
the shell does not oscillate radially but just expands or contracts exponentially. In both 
cases, the static shell is unstable up to the linear order. 

When ujr = 0, whether the AT-shell expands or contracts is determined by initial condi- 
tion. In this case, TZ\ can be expressed as 

TZ X = e WJ *-(7i. (4-75) 

Therefore, if Hi is initially positive, then C\ is positive and we find that the AT-shell expands 
exponentially. By contrast, if Tl\ is initially negative, then C\ is negative and it contracts 
exponentially. This behavior of the AT-shell can be understood as infinite period limit 
(T = 27t/wr — > oo) of the oscillation behavior f!4-74p . 

The appearance of the pure imaginary modes depends on the parameter y. By using the 
speed of orbital motion of each constituent particle measured by an observer who rides on 
the static AT shell, say v, the parameter y can be expressed as® 

y = r^- (4,76) 

It is easy to see that y is a monotonically increasing function of v. If the velocity v of each 
constituent particle is smaller than a critical value v c ~ 0.76 which corresponds to y ~ 1.4, 
the amplitude of the radial oscillation grows, whereas if v > v c , circumferential radius will 
expands or contracts exponentially with time. 

§5. summary and discussion 

We have investigated linear perturbation of the AT-shell around the static solution, and 
have found that the static state is unstable in the sense that, if the speed of orbital motion of 
the constituent particles (denoted by v ) is small, the perturbation of the shell's circumference 
radius oscillates with exponentially growing amplitude and the shell does not settle down 
into the static configuration, and if v is larger than the critical value v c ~ 0.76, the shell just 
expands or contracts exponentially with time. Whether the shell expands or contracts for 
large velocity depends on the initial condition. If initially the circumference radius is larger 
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than that of the static shell, then the AT-shell just expands exponentially. By contrast, if 
initially it is smaller than that of the static shell, the AT-shell just contracts exponentially. 

Since we invoked the numerical method to find the solutions, we have investigated a 
limited number of points in the domain < y < oo, not all y. But, as far as we have 
investigated, there exists at least one unstable mode for each y. Therefore, it is reasonable 
to conclude that the static AT-shell is unstable for all y up to 0(e). 

This result is compatible with the previous work given in Ref. [8] which showed that 
there exist momentarily static and radiation free initial states of the AT-shell which do not 
settle down into the static state, unless gravitational radiation extracts an infinite amount of 
C-energy from the AT-shell. More concretely, in this study, it was shown that, if the initial 
circumference radius of the AT-shell is greater than that of the expected final static state 
(having fixed a, A, R and ip s ), then the AT-shell cannot settle down into the static state, and 
if the initial circumference radius of the AT-shell is smaller than that of the static state, then 
it is not forbidden that the AT-shell settles down into the equilibrium static configuration. 
However, the existence of the unstable modes in our perturbation analysis does not depend 
on whether the initial circumference radius is greater than that of the static state or not. It 
means that, even if the initial circumference is smaller than that of the expected static shell, 
the AT-shell does not settle down into the static state. 

Hamity, Cecere and Barracd-® discussed the stability of the AT-shell and concluded 
that for y < 0.78049 (-Jy < 0.8836), the static configuration is stable and for y > 0.78049 
(yfy > 0.8836), it is unstable. Their result seems to be inconsistent to Ref. [8], and further, in 
the present analysis, we do not find any evidence for it. They investigated only the sequence 
of the static configuration by a Gedanken experiment. However, it should be noted that 
in order to obtain a definite conclusion on the stability of the static AT-shell, it would be 
necessary to solve the Einstein equations. In contrast to Hamity et al., we have solved the 
Einstein equations, although we have used linear approximation. 

Our perturbation analysis shows that the static AT-shell solution is unstable but does 
not say anything about the shell's finial state. To make it clear, perturbation theory will not 
be adequate and we will have to solve full equations by using numerical relativity. These 
will be future works. 
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